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ABSTRACT 

We present results for the QUEST RR Lyrae Survey at low galactic latitude, conducted entirely 
with observations obtained with the QUEST mosaic camera and the 1.0/1. 5m Jiirgen Stock 
Schmidt telescope at the National Observatory of Venezuela. The survey spans an area of 
476 deg 2 on the sky, with multi-epoch observations in the V, R, and / photometric bands 
for 6.5 x 10 6 stars in the galactic latitude range —30° ^ b ^ +25°, in a direction close to 
the Galactic Anticenter 190° ^ / ^ 230°. The variability survey has a typical number of 30 
observations per object in V and / and ~ 25 in R, with up to ~ 120 — 150 epochs in V and / and 
up to ~ 100 in R in the best sampled regions. The completeness magnitudes of the survey are 
V = R = 18.5 mag, and / = 18.0 mag. We identified 211 RR Lyrae stars, 160 bona fide stars 
of type ab and 5 1 candidates of type c, ours being the first deep RR Lyrae survey conducted 
at low galactic latitude. The completeness of the RR Lyrae survey was estimated in > 95 per 
cent and ~ 85 per cent for RRab and RRc stars respectively. Photometric metallicities were 
computed based on the light curves and individual extinctions calculated from minimum light 
colours for each RRab star. Distances were obtained with typical errors ~ 7 per cent. The RR 
Lyrae survey simultaneously spans a large range of heliocentric distances 0.5 ^ R ne i (kpc) ^ 
40 and heights above the plane — 15 ^ z(kpc) ^ +20, with well known completeness across 
the survey area, making it an ideal set for studying the structure of the Galactic thick disk. 

Key words: stars: variables: other, Galaxy: stellar content, Galaxy: structure, Astronomical 
Data bases: surveys 



1 INTRODUCTION 

In the Milky Way (MW), the thick disk hosts a very old (> 10 
Gyr) and relatively metal-poor ([Fe/H ] ~ —0.7) stellar population 
dWvsel2009l : lReddv & Lambertl2008l) . com prising ~ 10 per c ent or 
even up to ~ 20 of the thin disk mass (e.g. lJuric et ai]|2008l) . The 
Galactic thick disk therefore constitutes an important fossil record 
of the earliest stages of the formation of the Galactic disk, which 
could help understand the relevant mechanisms contributing to the 
formation of our Galaxy. From the perspective of galaxy formation, 
thick disks are remarkably relevant since these have proven to be an 
ubiquitous component of disk galaxies, as external galaxy surveys 
have shown that approximately 95 per ce nt of disk galaxies contain 
thick disks dYoachim & Dalcantorjbood) . 

On the other hand, despite many advances in the last few 
years, there is no consensus about the formation mechanism of 



cmateu@cida.ve 



the Galactic thick disk, or even on some fundamental proper- 
ties of its structure, mainly due to observational difficulties. Scale 
length measure ments for the MW thick disk range from Iir ~ 
2 to 4.7 kpc fcarolloe7aI1l2010l: [Larsen & Humphreys! 120031 : 
ICabrera-Lavers et alj l2005r . IChiba & Beeri |200C|) . which makes 
it unclear whether the thick disk is more radially extended than 
the thin disk or not, or how the density profile behaves at large 
radii. Another example are metal abundances, the iron abundance 
of thick disk star s seems to be fairly homogeneous in the verti- 
cal direction (e.g iKatz et alj|201 lLISoubiran et alj|2003f) . although 
IChen et alj feOllh recently report an appreciable metallicity gra- 
dient. Radial metallicity gradients have not been probed yet and 
there is also ongoing debate regarding the reach in [ Fe/H] at both 
the low and high metallicity tails of the thick di sk dBensbv et alj 



l2007l:lReddv & Lamberd|2008l : iReddv et alj|2005t) . In a recent pa- 
per. lBovv et al.l d201 lal lbl) argue against the need of a thin/thick disk 
decomposition, proposing that the Galactic disk can be described 
as a series of simple stellar populations having scale lengths and 
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Figure 1. Aitoff projection map showing the footprints of recent large-scale RRLS surveys in Galactic coordinates. Left: Deep surveys (V > 16), (blue) 
Present survey, (cyan) QUEST halo RRLS survey, (orange) SDSS Stripe 82, (green) LONEOS,(W) SEKBO, (purple) MACHO, (darkyellow) OGLE. Right: 
Shallow surveys (V < 16), (magenta) NSVS, (yellow) ASAS-3 and (light green) HATNET. The solid black line corresponds to the Celestial Equator. Survey 
characteristics are summarized in Table[T] 



scale heights which vary smoothly as a function of metallicity and 
a— element abundance. This illustrates how the issue of Galactic 
thick disk structure is far from settled and reliable measurements 
on structural, kinematic and chemical properties are necessary to 
gain a better understanding of the MW's thick disk which in turn, 
will impose strong constraints on the possible contribution of the 
different formation mechanisms. 

Observationally, disentangling the structure of the Galactic 
disk is a challenging task since the Sun is embedded in it. Numer- 
ous studies have been conducted at high galactic latitudes (\b\ > 
30°) in order to study stars at a height above the Galactic plane 
z ~ 1 — 3 kpc such that the thin di sk density has decayed suffi- 
ciently, reducing contamination (e.g.|Oiha 2001; Siegel et alj2002t 
iBrown et al.ll2008l ; IWvsell2009t) . but also limiting the coverage in 
the radial direction. At intermediate to low latitudes (\b\ < 30°) the 
problems of crowding and the high and variable extinction become 
important, and the contribution of the thin disk has to be modelled 
dRobin et alJl996|lLarsen & Humphrevsl2003l ; ICarollo et alj201Ct 
IdeJong et alkoiOl) . 

RR Lyrae stars (RRLS) have several fundamental properties 
that offer advantages in dealing with these problems inherent to low 
galactic latitude observations. RRLSs are luminosity and colour 
standards. Since RRLSs are horizontal branch stars, they have a 
well known absolute magnitude which shows little spread jSmithl 
1 19951) . and allows the computation of di stances wi th sm all uncer- 
tainties (< 10 per cent. lVivas et alj|2004 hereafter fvoj) . Adition- 
ally, during the phase of minimum light of the pulsation cycle, 
RRLSs of type ab have approximately the same effe ctive temper- 
ature, and thus show very little dispersion in colour dSturchll 19661 ; 
iDav et al. 1 120021) . This property enables the use of RRab stars as 
colour standards to measure extinctions up to the distance of each 
individual star; which is a crucial point in a low latitude survey 
since reddenning changes drastically along different lines of sight 
and has a strong dependence with distance. Also, since RRLSs 
trace old (>10 Gyr) and mainly metal-poor ([Fe/Hj < —0.5) stel- 
lar populations |Smith 19951; iDemarque et alj|200fj|). the expected 
conta mination from the thin disk is negligible dMartin & Morrison! 
Il998l) . Finally, RRLSs are pulsating stars with relatively short pe- 



riods (0.3 — 1.0 d) and large light curve amplitudes (V ~ 0.3 — 1.2 
mag), which makes th em easily identifiable by a pho tometric multi- 
epoch survey (see e.g. lV04lKinemuchi et alj|2006h . 

RRLSs have been extensively used as tracers of the Galactic 
halo by numerous surveys which have studied their spatial distri- 
bution in a wide range of distances, from very near the Galactic 
centre up to large distances R^ a [ > 100 kpc. This is illustrated in 
Figure[T]which shows, in an Aitoff projection map in galactic coor- 
dinates, the footprint of recent large-scale RRLS surveys including 
the present one. Table[T]summarizes the characteristics of each sur- 
vey. 

Unlike in the Galactic halo, the distribution of RRLSs in 
the thick disk has been less thouroughly studied, in particular 
due to the lack of deep RRLS surveys at low/intermediate galac- 
tic latitudes. As illustrated in Figure Q] the Galactic Plane area 
\b\ < 20° has only been covered b y the large-scale yet shallow 
(V < 14 or R h „i < 5 kp c) ASAS-3 dPoimanskill2002h and NSVS 
dKinemuchi et al .1 120061) surve ys, as well as t he co mpila tions of 
nearby RRLS (V< 13) from ILavdenl dl994 1 19951) and iMaintzl 
d2005l) ; having been avoided by deep surveys (16 < V < 20). 
This makes ours, the first deep large-scale RRLS survey con- 
ducted at low Galactic latitudes. Furthermore, our survey probes the 
out er regions of the thick disk, outside the s olar circle. In particu- 
lar.lLavdenl dl995h ; lMaintz & de Boen d2005l) and lKinemuchi et al.l 
d2006T NSVS) have studied the distribution of thick disk RRLSs in 
the vertical direction, but the limited coverage at low galactic lati- 
tude prevented the determination of the scale length as well as the 
exploration of deviations from an exponential profile such as the 
warp and flare, which have been ob served in the Galactic thin and 
thick disks with other tracers (e.g. Lopez-Corredoira et alj|2002t 
iMomanv et al ] |2006l ; lHammerslev& L6pez-Corredoirall201l|) . For 
these reasons it was essential to conduct a deep RRLS survey at low 
galactic latitude, which provided an ample spatial coverage both in 
the radial and vertical directions. 

In this paper we will present the catalogue of the QUESlO] 
RR Lyrae survey of the thick disk, which spans a total area of 476 
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Table 1. Characteristics of recent large-scale RRLS surveys 



Survey Filters Area (deg 2 ) Completeness Telescope Observatory Reference 



QUEST 


VRI 


476 


V < 19 


1.0 m Schmidt 


NOV 


This Work 


QUEST (V04) 


V 


380 


V < 19 


1.0 m Schmidt 


NOV 


V04 


NSVS 


ROTSE-NT 


~ 31000 


V < 14 


4 x 200 mm ROTSE-I 


Los Alamos 


Kinemuchi et al. (20061 


ASAS-3 


V 


~ 31000 


V < 14 


2 x 200 mm ASAS 


Las Campanas 


Poimanski ("20021 


LONEOS 


LONEOS-NT 


1430 


V < 18 


0.6 m Schmidt 


Lowell 


Miceli et al. (20081 


SEKBO 


BmRm 


1675 


V < 19.5 


1.27 m MACHO 


Mount Stromlo 


Keller et al. (20081 


HATNET 


I 


67 


V < 15 


1 1 cm HATNET 


Fred L. Whipple 


Hartman et al. (20041 


OGLE 


1 


0.87 


V < 20.5 


1.3 m Warsaw 


Las Campanas 


Udalski et al. (1998) 


SDSS Str82 


ugriz 


249 


V <22 


2.5 m SDSS 


Apache Point 


Watkins et al. (20091 














Sesar et al. (20101 


MACHO 


BmRm 


~75 


V <20 


1.27 m MACHO 


Mount Stromlo 


Kunder & Chabover (20081. 














Alcock et al. (20031 



deg 2 at -30° < b < +25°, with multi-epoch V,R,I observations 
obtained with the 1.0/ 1.5m Jiirgen Stock Schmidt telescope and 
the QUEST mosaic camera, at the National Astronomical Observa- 
tory of Venezuela (NOV). In Section[2]we describe the observations 
and spatial and time coverage of the survey, as well as the photo- 
metric processing, error estimation and calibration. In Section[3]we 
describe the procedures used in the identification of RRLSs, we 
estimate the completeness and possible contamination of the sur- 
vey and we describe the computation of periods, amplitudes, red- 
denings, photometric metallicities and distance for the identified 
RRLS. In Section [4] we summarize the characteristics of our sur- 
vey. The computation and analysis of RRLS density profiles for the 
Galactic thick disk and halo will be presented in an upcoming paper 
of the series. 



2 THE SURVEY 
2.1 The Data 

The present survey spans a total area of 476 deg 2 and makes use 
of archive observations obtained between late 1998 and mid-2008 
with the QUEST camera and the 1.0/1.5 m Jiirgen Stock Schmidt 
telescope at the NOV in Llano del Hato, Venezuela. The camera 
is a 16 CCD mosaic, with 2048 x 2048 pixel CCDs layed out in a 
4x4 array, having a field of view of 2P3 x 2°3 and a scale of If 02 
pixel -1 . The QUEST camera was designed to operate optimally 
in drift-scanning mode and can be fitted with up to four different 
filters at a time. For a more det ailed descript i on of the c amer a and 
the drift-scanning technique see iBaltav et alj J2002h andlvol. 

For the present survey we used archive drift-scan observations 
at low Galactic latitudes Qb\ < 30°) approximately in the direc- 
tion towards the Galactic anticentre. These observations were ob- 
tained by different observational projects with very different goals, 
but since the equipment and method of o bserva tions are the same 
as the original QUEST RR Lyrae survey Jvoli , we have decided 
to keep that name for the present survey. Many of the observa- 
tions come indeed from the QUEST (quasar) survey (Baltay et al. 
2002; Rengstorf et al. 2004, V04), and there is also an important 
number of scans from the CID/Q Equatorial Variability Survey 
tericeno et a"l]|2003h . The latest spans the equatorial region in the 
full range of right ascension and has thoroughly surveyed the Orion 
star forming region, resulting in a large number of observations in 
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this area l lBriceno et alj|2005l iDownes et alj|2012l ). Our sample is 
restricted to observations in the range 60° ^ a ^ 135° in right as- 
cension (corresponding to the range 190° ^ I ^ 230° in Galactic 
longitude), where the Galactic Plane crosses the celestial equator 
and drift-scan observations in the declination range |8| ^ 6.5° are 
available. 

The observations were obtained using different filter combi- 
nations and have different spatial coverage and cadence, all of 
which changed over time, resulting in a highly inhomogeneous 
sample. The survey is composed by 452 drift-scans correspond- 
ing to 19,238 hours of observation and 6.3 Tb of data in the V, R 
and / photometric bands, obtained in a total of 427 nights between 
December 1998 and June 2008. Although ideally an RRLS survey 
would use a blue filter, since the amplitudes of variation of RRLS 
would be larger, we did not include these observations in our sam- 
ple since these were very few due to the low sensitivity of QUEST 
CCDs in this band. 

The present survey include s the low galactic latitude observa- 
tions from the RRLS survey of IV04l which consisted of QUEST 
drift-scans at 8 = — 1° obtained between December 1998 and April 
2001. This results in an overlap area of 90 deg 2 at both ends of 
the present survey (see Figure QJ in the right ascension ranges 
60° sC a s: 90° and 120° < a < 135° and the declination range 
—2° ^8^0°. New observations (post-2001) were added in this 
area and the existing ones from lV04| were reprocessed using the 
procedures described in the following section, which were devised 
to make an optimal use of data in the three available filters V,R,I. 
Finally, data in the ran ge 60° ^ Ot ^ 90° ar e also part of the CIDA 
Deep Survey of Orion jDownes et al JI2012L in preparation). 

2.2 Data Reduction and Photometry 

In the following sections we will describe the procedures used in 
the object detection, astrometric reduction, aperture and PSF pho- 
tometry. 

2.2. 1 Object Detection, Astrometry and Aperture Photometry 

Each individual drift-scan observation was reduced using Offline, 
the standard pipe line developed by the QUEST colaboration 
( IBaltav et alj|2002f) . The Offline pipeline automatically performs 
basic reductions, object detection, aperture photometry and astro- 
metric reductions for each image frame in a drift-scan observation. 
This processing was conducted independently for the images ob- 
tained by each of the camera's CCDs , resulting in 16 independent 
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catalogues having position in celestial equatorial coordinates, pho- 
tometric magnitude and error for each point-source, as well as pho- 
tometric parameters such as FWHM, ellipticity, sky background 
flux and error, etc. 

Image reduction consists in standard bias, dark and flat- 
fielding corrections. Object detection is made at the 3-<J level and 
aperture photometry is performed over all detected point-sources 
with an aperture radius of one FW HM. The ast rometric matrices 
are computed using the USNO-A2 (Monet 1998) astrometric cata- 
logue as a reference. The astrometric solutions computed by Offline 
have a typical accuracy of ~ 0"15, which is sufficient for the cross- 
identification of an object in the different catalogues. 



2.2.2 PSF Photometry 

Since our drift-scan observations include very crowded regions at 
low Galactic latitude, the use of PSF photometry was necessary. 
Figure[2]shows a plot of the typical number of stars per image in the 
V, R and / bands, as a function of right ascension a for a drift-scan 
at 5 = —3° which spans the entire right ascension range covered by 
this survey. The right-hand vertical axis shows the scale in terms of 
Nc the number of stars in a circle with a 16" radius (which is twice 
the aperture diameter in an observation with a bad seeing of 4"). 
The plot illustrates how the number of objects per image increases 
when approaching the Galactic Plane. In particular Nc > 1 in the 
range 95° < a < 1 15°, meaning that the typical separation between 
two objects in an image is less than twice the aperture diameter and, 
therefore, aperture photometry for most stars will be contaminated 
by neighbours. Therefore, the computation of PSF photometry was 
restricted to the right ascension range 90° ^ a ^ 120°, in which 
Nc > 0.8 (indicated in Figure [2] with dashed black lines). In low 
density areas, consistency was checked by computing both aperture 
and PSF photometry. 

PSF phot ometry was obtained using the standard DAOPHOT 
dStetson|[l987h routines implemented in tied together in 

an automated pipeline called OfflinePSF, which is customized for 
QUEST drift-scan observations. 

OfflinePSF includes as a first step, prior to the computation 
of PSF photometry, the recentering and recomputation of aperture 
photometry for all point-sources previously detected with Offline. 
This is necessary since initial Offline detections have centering er- 
rors of ~ 1 pixel, which are reduced to ~ 0. 1 pixel by using PHOT's 
recentering options. Using the recalculated centroids and aper- 
ture photometry magnitudes, PSF photometry is obtained accord- 
ing to an iterative scheme where, for each individual image, input 
parameters such as the fit radius r« ; = 1 x FWHM and the model 
PSF radius r ps f = 3 x rjj, are computed independently. Then, from 
90 to 120 PSF stars are selected, uniformly distributed accross the 
image, without neighbours within a radius r ps f and, for distances 
less than 2r„ s f requiring neighbourghs be at least one magnitude 
fainter and there are no bad pixel columns near. An initial PSF 
model is made with the DAOPHOT PSF task, with a second or- 
der spatial variation across the image. The PSF model is fit around 
PSF stars with the NSTAR and SUBSTAR tasks, PSF stars with 
newly detected neighbourghs are discarded and the PSF model re- 
computed. The PSF model is then fit to all point-sources in the en- 
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Figure 2. Number of stars per image as a function of a, for a drift-scan 
centered at S = —3°, with filters V, R e / (*, x , + respectively). The right- 
hand vertical axis indicates the typical number of stars Nc inside a circle of 
radius 16" (C(16")). Dotted grey lines indicate the regions corresponding 
to galactic latitudes b = ±20° , ±40° . The black dashed lines delimit the 
region were PSF photometry was computed. 



tire image using ALLSTAR and new objects are detected if present 
in the resulting subtracted frame, iterating twice more over these 
two steps. 

The Offline and OfflinePDF processes result, for each drift- 
scan, in a photometric catalogue per CCD, containing the equatorial 
coordinates a, 8, instrumental aperture and PSF magnitudes and er- 
rors, for each point-source object detected in the full drift-scan. 

2.2.3 Photometric Calibration 

The photometric calibration was done in two stages: the calibration 
of a master catalogue by means of second order photometric stan- 
dards and the normalization of individual drift-scan observations to 
this master photometric catalogue. 

The master catalogue is composed of 3 to 4 individual scans 
per declination stripe, observed with good and constant seeing and 
sky background throughout the whole duration of the scan. The 
combination of these scans allowed for a much more homogeneous 
master catalogue, than would have been obtained if only one scan 
had been used, due to empty regions caused by bad columns and 
defective regions of the CCDs. The master catalogue obtained con- 
tains 6,513,705 objects and its coverage is illustrated in the density 
map shown in Figure[3] 

The instrumental magnitudes for each object in the master cat- 
alogue were calibrated using second order photometric standards. 
A total of 4526 photometric standards were used, spanning the 
magnitude range 10 < V < 17 with colours —0.2 < V—R < 1.2 
and identified as non-variable stars in any one of the VRI bands, 
with a spatial distribution illustrated in Figure[3] Photometric stan- 
dard stars ind i cated by black dots (2621 in total) were set up by 
iDownes et al. I d2012h using observations obtained with the 1.2 m 
telescope and the Keplercam camera, at the Fred L. Whipple Ob- 
servatory, Mt. Hopkins, US. The asterisks indicate the_28 stan- 
dards having 8 = -1° and a < 175°, established by |V04| with 
observations from the 1.0m YALC0 telescope at the Cerro Tololo 
Interamerican Observatory (CTIO), Chile. Finally, as part of this 
work, we established 1877 second order standards in the range 
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Figure 3. Density map of stars in the master photometric catalogue and 
spatial distribution of second order photometric stantards used for photo- 
metric calibration. The gray scale represents stellar density as indicated by 
the colour bar, per 0.24 deg 2 bin (equivalent to half the area of a QUEST 
CCD). The dots indicate the positions of sec ond order photometri c stan- 
dards established by different authors: (blue) iDownes et alj 120121) . (red) 
|V04| and (orange) this work. 



—4° ^ 8 ^ —2° (denoted in Figure [3] with dark gray dots), us- 
ing observations obtained during December 2004, with the 0.9 
m SMARTS^ telescope at CTIO. Data for these 1877 second or- 
der standards is summarized in Table [2] The photometric calibra- 
tion was performed independently in declination stripes having the 
width of a QUEST CCD and in 1 -hour-long intervals in right as- 
cension where the distribution of standards allowed it. The cata- 
logue regions centered around 8 = — 1° with 110° ^ a ^ 135° and 
8 < — 5° with a > 130° lacked photometric standard stars, so these 
were calibrated using l ong re ference scans which ext ended as far as 
the area covered by the(V04j standards (asterisks) and lDownes et al.1 
(black dots) respectively. 

The photometric solutions for each scan were obtained using 
zero-point and V — R or V—/ colour terms, depending on the filters 
available for each observation. A mean photometric precision of 
~ 0.021 mag was achieved for the V,R and / bands. 

2.2.4 Normalization 

The master photometric catalogue, described in the previous sub- 
section, provides the calibrated reference photometry for all survey 
objects. The individual observations, used to construct time series 
for each object (see Sec. |3,2t used in the RRLS search, were cali- 
brated by normalizing each scan to the mast er ref erence catalogue. 

The normalization process, taken from lV04l consists in iden- 
tifying all objects of a particular scan in the reference master 
catalogue and computing the mean magnitude residuals in both, 
throughout the whole scan in 0'25 bins in right ascension. Since 
most stars (> 95 per cent) are non-variable and the number of stars 
in each bin is sufficiently large (~ 1500 — 3000), the average resid- 
uals obtained result in a precise estimation of the zero point needed 
to calibrate each bin of the scan. This normalization process allows 
for the correction of systematic trends in the magnitude residuals 
caused by changes in the atmospheric conditions during each ob- 
servation, which in our case is critical given that a single drift-scan 
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Figure 4. Magnitude standard deviation a as a function of mean V mag- 
nitude, for stars in a typical 0.5° x 15° region centered around (a, 8) = 
(85°, +3°), with a mean number of ~ 25 observations per star. The solid 
line indicates the polynomial fit of the average o" in 0.5 mag bins, as a func- 
tion of mean magnitude. Error bars indicate the standard deviation around 
each of these averages. 



can be up to 7 to 8 hours long (see Fig ure 3 in lV04l for an example 
of the normalization process). The typical resulting standard devia- 
tion of the residuals in each bin was ~ 0.015 mag, giving zero point 
errors ~ 0.002 mag. 

2.2.5 Photometric Errors 

Taking advantage of the fact that each star in our survey has mul- 
tiple observations in each filter, it is suitable to use the standard 
deviation a of magnitudes of non-variable stars as the representa- 
tive photometric error in the given filter. Since the majority of stars 
in any given field are non-variable, these give rise to the main trend 
in a plot of a versus mean magnitude, as seen in Figure [4] 

In the figure the crosses show the average standard deviation 
computed as a function of mean magnitude in 0.5 mag bins using an 
iterative 3o— clipping. The sharp increase in o" for objects brighter 
than V ~ 14.5 mag is due to objects saturated in most observations. 
This behaviour provided us with a means to empirically determine 
the saturation magnitude and, for each object, remove the individ- 
ual observations brighter than this limit. The solid line indicates a 
polynomial fit of the average a versus magnitude (crosses) and the 
error bars correspond to the standard deviation around this aver- 
age. The polynomial fit consequently provides the typical magni- 
tude standard deviation due to photometric uncertainties only, as a 
function of an object's mean magnitude, which we assumed as its 
photometric error. 

For each survey star we computed the corresponding photo- 
metric error as the the typical magnitude standard deviation given 
by the polynomial fit, for the corresponding mean magnitude. This 
error analysis was conducted independently for observations in 
each of the V,R and / filters and dividing the entire survey in sepa- 
rate regions of A8 = 0.5° times Aa = 15°. 

2.3 The Photometric Catalogue 

The final photometric catalogue spans a total area of 476 sq. deg. 
having 6,513,705 objects with multi-epoch observations in the V, 
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Table 2. Magnitudes for second order photometric standards. 



ID a 8 V R I 





(°) 


(°) 


(mag) 


(mag) 


(mag) 


ss4()00 
ss4001 
ss4002 
ss4003 
ss4004 


82.401749 
82.420319 
82.431961 
82.436493 
82.437607 


-2.04192 
-2.10225 
-2.04361 
-2.11058 
-2.03612 


13.446 ± 0.003 
14.386 ± 0.005 
16.905 ± 0.020 
16.630 ± 0.018 
16.784 ± 0.019 


12.238 ± 0.003 
13.201 ±0.005 
14.526 ± 0.009 
15.197 ±0.014 
15.329 ±0.015 


12.822 ± 0.003 
13.800 ± 0.004 
15.849 ± 0.012 
15.912 ±0.012 
16.046 ±0.013 



Note: Table|2]is published in its entirety in the electronic edition of the journal. A portion is shown here 
for guidance regarding its form and content. 
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Figure 5. Spatial coverage of the survey in the V (top), R (middle) e / (bottom) bands. The colourbars indicate the mean number of epochs per object in each 
filter. 



R, and / bands. Table [3] summarizes the saturation (m s ), complete- 
ness (m c ) and limiting (m/) magnitudes in the three photometric 
bands, as well as the typical errors corresponding to these magni- 
tudes. 

The spatial coverage of the survey as well as the typical num- 
ber of epochs per object, in each of the photometric bands, are il- 
lustrated in Figure[5]in equatorial coordinates. As can be seen in the 
figure, the number of observations per object goes from ~ 10 up to 
~ 150 in the V and / bands and up to ~ 60 in the R band. The num- 
ber of epochs per star per filter is quite inhomogeneous through- 
out the survey, as illustrated by Figure[5] nevertheless most objects 
have at least 10 observations in two photometric bands; while ob- 
jects in the thoroughly sampled Orion region have more than ~ 30 
observations in each of the V, R and / bands. 



3 THE RR LYRAE SEARCH 
3.1 Light Curve Simulations 

The present survey has very inhomogenous time sampling and 
number of epochs in the different filters, which in turn change in 



Table 3. Saturation, completeness and limiting magnitudes of the sur- 
vey. 



Filter 


m s 


Am, 


m c 


Am ( . 


111/ 


Am/ 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


V 


14.0 


0.018 


18.5 


0.064 


19.70 


0.15 


R 


14.0 


0.018 


18.5 


0.072 


19.75 


0.15 


I 


13.5 


0.020 


18.0 


0.078 


18.80 


0.15 



different regions of the survey. These differences can have a size- 
able impact in the completeness of the resulting RRLS sample, 
which must be properly characterized. We used light curve sim- 
ulations as a tool to test and fine-tune the variable star detection 
and period searching techniques used to identify RRLS stars (Sec. 
13.21 and \3. 3. Il l, as well as to characterize the completeness of the 
RRab and RRc samples as a function of position in the survey 
(Sec. l3"33t . 

A total of 425,000 synthetic light curves were generated for 
RRLS of each type, ab and c, and the same number for non- 
variable stars. The time sampling and photometric errors as a func- 
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Table 4. Parameters used in the simulation of RRflfo and RRc light curves. 



Parameter 


RRab 


RRc 


Reference 


(P) 


0.539 d 


0.335 d 


V04 


ol> 


0.09 d 


0.07 d 


V04 


(AmpV) 


1.04 mag 


0.536 mag 


V04 


OAmpV 


0.24 mag 


0.13 mag 


V04 


AmpS 


0.74064 x AmpV +0.0361 


Alcock et al. (2003) 


Amp/ 


0.66800 x AmpV +0.0501 


Dorfi & Feuchtinger (1999) 


(V) 


13.0- 


-20.0 




((V-R),a v - R ) 


(0.30;0.26) 


Alcock et al. (2003) 


«V-/),a V -7) 


(0.42; 0.20) 


Alcock et al. (2003) 



tion of magnitude (see Sec. 12.2.51 were taken from actual objects 
uniformly distributed throughtout the entire survey area. This en- 
sures a truly representative realization of the time spacing, the num- 
ber of observations per filter and the number of quasi-simultaneous 
observations in different filters, for each synthetic object. This point 
is of particular importance given the differences both in the cadence 
and spatial coverage in the different photometric bands as a func- 
tion of position, as illustrated in Figure [5] 

The parameters used to create synthetic light curves for 
RRa£> and RRc stars are summarized in Table|4] We generated syn- 
thetic stars with random mean V magnitudes (V) uniformly dis- 
tributed in the range [13.0,20.0]. For each object, the correspond- 
ing mean magnitudes (R) and (/) were computed from the mean 
V — R and V — I colours, randomly drawn from gaus sian distribu- 
tions w ith means and standard deviations taken from I Alcock et al.l 
d2003l) . indicated in Table|4] The pulsation period P of each object 
was randomly drawn from a gaussian distribution and the K-band 
amplitude computed from the (O osterhoff I) Period- Amplitude re- 
lation from lCacciari et al" I d2005h . with random gaussian noise with 
sigma 0.12 mag, for RRab stars. For RRc stars, since the Period- 
Amplitude relation is weak, the V-band amplitude was randomly 
drawn from a gaussi an dis tribution assuming a mean a standard 
deviation taken from |V04 see Table [4] This way our simulated 
population of RRLS shows the period-amplitude trend for Oost- 
erhoff I RRab stars ( which represent ~ 80% of all RRafc stars in 
the Galactic halo, e.g. ISesar et al] |2010j), albeit not the Oosterhoff 
dicotomy. Note that although we do not simulate the Oosterhoof 
dicotomy, nor the Blazhko ef fect (an ampli tude modulation known 
to occurr in some RRLS, e.g. Kovacs 2009), we do not expect they 
will affect our completeness estimates for the survey. For a given 
period, Oosterhoff II stars would have even larger amplitudes than 
Oosterhoff I ones making them even easier to identify. On the other 
hand, the amplitude of RRLS is large enough such that the (rela- 
tively small) variations due to Blazhko effect will only translate in 
the light curves appearing slightly more noisy. 

Given the V-band amplitude (AmpV), the corresponding am- 
plitudes in the R and / bands were computed acc ordi ng to the lin- 
ear rel at ionships from [jporfi & F euchtingeil l l 19991) and lAlcock et all 
d 19991) ; I Alcock et alj dlooll) respectively, including the gaussian 
uncertainty estimated for these relationships. This procedure cor- 
rectly reproduces the ampli tude ratios observed for RRLS in differ - 
ent photometric bands (e.g. lSmith : 1995; Sterken & Jaschek] |2005l) . 
Finally, an initial phase offset randomly drawn from a uniform dis- 
tribution in the G [0, 1) range, was added to the synthetic light 
curve of each object which in t urn was s ample d fro m RRafc or 
RRc light curve templates from iLavdenl d 19981) and Ivivas etafl 



d2008l) respectively, in order to obtain the magnitudes correspond- 
ing to each epoch of observation (in each filter). A random gaus- 
sian photometric uncertainty was added to each of the synthetically 
produced magnitudes, using the error vs magnitude curves for each 
filter derived for the survey (see Sec. 12.2.51 . 

Time series data were also generated for synthetic non- 
variable stars with (V) € [13.0, 20.0]. Mean (R) and (/) magnitudes 
were computed from V — R and V — I colours, randomly drawn 
from uniform distributions in the ranges [—0. 1,1.5] and [—0.2, 3.0] 
respectively, which correspond to the entire V — R and V —I colour 
range for main sequence stars according to iKenvon & Hartmannl 
dl995h . 



3.2 Variable star identification 

The first step in the search for RRLS consisted in identifyi ng vari- 
able s tars, by means of the variability indices defined by IStetsonl 
dl996h . 

The J and K vari ability indices were initially defined by 
Welch & Stetson] { 993b and later modified and combined by 
Stetson! dl996T) n the definition of the L index, specifically designed 
for the automated identification of pulsating variables and the re- 
duction of contamination due to spurious variables caused by de- 
viant photometric observations. The definition of the L index uses 
the fact that, for pulsating stars, quasi-simultaneous observations 
in two photometric bands must experience deviations (with respect 
to the corresponding means) which show a positive correlation, as 
opposed to non-variable stars for which this deviations would be 
uncorrelated. 

The QUEST drift-scans comprising this survey have data in 
up to four different filters (see Sec. [23} obtained almost simultane- 
ously, since the mean separation between observations in two dif- 
ferent filters is At ~ 2 — 6 min (given by the time it takes an object 
to cross the different CCD rows of the camera), which compared 
to the typical period of an RRLS star (P ~ 0.55 d) can be consid- 
ered negligible {At IP ~7 x 10~ 3 ). This makes it an ideal set of 
observations to use lStetsonl 's variability indices. 

Ca ndidate variable stars were identified based on the three 
IStetsonl L indices that can be computed with our observations, i.e. 
LyR,Lyj and Lrj. Given that the survey coverage and sampling are 
very inhomogeneous in the different filters (see Figure[5). the crite- 
ria used in identifying variable stars needed to be sufficiently flex- 
ibe in order to consider e.g. objects which have observations in any 
two bands and not the third, yet without seriously reducing the com- 
pleteness in the process. Therefore, we selected as candidate vari- 
able stars, those objects having at least one of the three possible L 
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Figure 6. Left: Lvr versus V magnitude, for synthetic non-variable (black), RRab (*) and RRc (triangle) stars. The dashed line indicates the threshold chosen 
for the selection of variable stars. Right: Completeness of the variable star identification, estimated from synthetic light curves of RRab (blue) and RRc (green) 
stars. The dotted lines indicate the survey's saturation and completeness magnitudes in the V-band (see TablefJ} 



indices greater than a thresh old value of 1, slightly higher than the 
threshold value 0.9 used bv lStetsorJ ( fl996h in the identification of 
Cepheid variables. This threshold value was chosen based on re- 
sults from light curve simulations of RRab, RRc and non- variable 
stars with the survey's time sampling. 

A plot of Lvr as a function of V is shown in Figure [6] (left). 
As expected, the larger amplitude RRab stars (blue dots) have on 
average large L values than the lower amplitude RRc stars (green 
dots), while non-variable (black dots) stars have L values around 0. 
The criteria used to select variable star candidates (at least one L 
index higher than 1, see Sec. I3.2| l ensures there will be negligible 
contamination from non-variable and low-amplitude (< 0.2 mag) 
variable stars, while still recovering most variable stars, regardless 
of whether these have a small number of observations (< 10 — 12) 
in some filters. The completeness of the variable candidate sample 
achieved using this criteria is illustrated as a function of mean V 
magnitude in Fig[6](right). Within the saturation and completeness 
magnitude limits (dashed lines), the completenetess in the variable 
star identification was estimated in 99. 1 per cent for RRab and 98.3 
per cent for RRc stars. 

Using this criteria we obtained a sample of 486,766 candidate 
variable stars, corresponding to 7 per cent of the total number of 
stars in the survey. 



3.3 RR Lyrae identification 

3.3.1 Period Search 

The search for RRLS was conducted using a version of the 
lLafler & Kinmanl i 19651) met hod for the ide ntification of periods 
of variable stars, modified bv lStetsonl j 19961) and extended in the 
present work in order to incorporate the information of multi-filter 
time series data. 

The method consists in computing, for a series of trial periods, 
the parameter 



S(M): 



lf M co(M+l)lm ; -m i+1 | 



(1) 



where m, and m,+i are consecutive magnitudes in the M filter, in 



the light curve folded with the trial period, Nm is the number of ob- 
servations and co(/, i+ 1) are weights proportional to the magnitude 
errors o,-, given by 



C0(i',i'+1) 



or + cx 



(2) 



;+l 



The true period corresponds to the minimum S(M), since the 
light curve is smooth when folded with the correct period and there- 
fore the sum of differences between consecutive magnitudes is min- 
imal. Since S(M) considers only the data in the M filter, in order to 
incorporate the information available for each object in all three 
V, R and / filters we used the comp osite index Svri as t he average 
of the 5(M) for each filter (see e.g. lWatkins et alj|2009l) . weighted 
by the corresponding number of observations, as given by 



Svri 



N V S(V)+N R S(R)+N,S(I) 
Nv+Nr+Ni 



(3) 



where Nv,Nr,Nj correspond to the number of observations in 
each of the V,R,I filters respectively. By minimizing the Svri pa- 
rameter given by Equation [3] we are using the data in all the avail- 
able filters, requiring that the best period be the one that smoothes 
the period-folded light curves in the three photometric bands. 

RRLS are known to exhibit a narrow range in temperatures, 
and therefore colour, which is why the use of colou r cuts is fre- 
quent in RRLS searches (e.g. llvezic et al feoollvol) . However, in 
low galatic latitude surveys, stars are severely affected by strong 
and highly variable extinction, and the use of colour cuts is better 
avoided in order to ensure the survey's completeness, albeit increas- 
ing the computational cost of conducting the period search over a 
much larger sample. Therefore, the period search was conducted 
over the full sample of 486,766 variable candidates selected with- 
out imposing any sort of colour cuts. The search spanned the period 
range from 0.2 — 1 .2 d with a 10~ 4 d period step and for each star, 
and a second order search was performed with increased period res- 
olution around the 5 best periods initially found, using a period step 
of 10~ 5 d. The best 3 periods found after the refined search were 
reported for each object characterized by a statistical significance 
parameter A defined as 



A= S VR! (P)-(S VRI ) (4) 

where (Svri) and o(Svri) are the mean and standard deviation 
of the distribution of Svri values for all trial periods. 

3.3.2 Period Aliases 

The emergence of alias or spurious periods is a common problem 
when trying to identify the period of a signal from time samplings 
which may or may not be uniform. In general, these can simply be 
harmonics of the true period or they can come up as a consequence 
of an underlying periodicity in the ti me sampling, being there fore 
external to the signal itself. Following lLafler & KinmarJ i ll 9651) we 
can express the relationship between the true period P and spurious 
periods n as 

where A; is a rational number, which produces the true period har- 
monics; and p is t he ex ternal period of the sampling. As noted by 
lLafler & Kinmanl d 19651) the external periods p of interest are those 
which produce period aliases in the expected range of the signal, in 
this case in the period range of RRLSs, since the rest can be ruled 
out on physical grounds. 

In the present case of the RRLS search, the typical one day 
periodicity between the drift-scan observations composing our sur- 
vey, leads to the emergence of spurious periods in the range 0.2 < 
P(d) < 1 expected for RRLSs. This spurious periods, commonly 
termed one-day aliases, are obtained from Equation[5]with p = ±1 
and k = 1 . Other period aliases can arise due to underlying period- 
icities in the time sampling, as a consequence of the inhomogeneity 
in the time sampling of our survey. 

In order to empiricaly identify the most frequent alias periods 
in our sample, we conducted the period search described in the pre- 
vious section, over the synthetic time series data for ~ 130,000 of 
our simulated RRab and RRc stars (see Sec. 13. II Figure [7J shows 
plots of the recovered period as a function of the true period for 
synthetic RRab (left) and RRc (right) stars. The most populated 
sequence in the plot is the identity, indicated by the red shaded 
area. Different alias periods given by Equation [5] match the other 
sequences observed in Figure [JJ The most prominent loci defined 
by the alias periods identified in the sample are shaded in differ- 
ent colours, as detailed in Table [5] The width of the shaded loci 
in Figure [7J indicates a separation smaller than 2 per cent from the 
corresponding main locus. 

Figure[7Jshows that the majority of the periods lie in the iden- 
tity locus. If the recovered period differs in less than 10 per cent 
from the true period, we consider it correctly recovered, with this 
definition we find that for RRab 80 per cent of the times the period 
is recovered, while for RRc stars it is only ~ 46 per cent. For the 
majority of the remaining stars the recovered period is an alias of 
the true period. This indicates, as can be seen in the figure, that the 
fraction of RRc stars affected by spurious periods is much larger 
than that for RRab stars. The frequency of occurence of the alias 
periods indicated in Figure [JJ is summarized in Table [5] for both 
types of RRLSs, computed from ~ 130.000 simulated RRab and 
RRc light curves. As can be seen both in the tabl e and in Figure [7J 
one-da y aliases are the most frequent as noted bv lLafler & Kinmari 
j 19651) . Nevertheless, other spurious periods can arise, albeit less 
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Table 5. Occurrence frequency of alias periods. 



Type 


k 


P 


Frequency (per cent) 


Color 






(d) 


RRab 


RRc 




1-d alias 


1 


+1 


1 1.9 


13.1 


Green 




1 


-1 


2S.4 


36.2 


Green 


1/2-d alias 


1 


+ 1/2 


1.9 


1 6 


Cyan 




1 


-1/2 




6.1 


Cyan 


1/3-d alias 


-1 


1/3 


12.6 




Orange 


1/4-d alias 


-1 


1/4 


17.4 




Grey 


Harmonics 


2 




8.4 


7.7 


Blue 




1/2 




1.7 


4 


Blue 


1-d alias 


1/2 


-1/2 




15.4 


Yellow 


harmonic 












1/2-d alias 


1/2 


-1/4 


17.7 


19.4 


Magenta 


harmonic 













frequently, such as the 1/3-day and 1/4-day aliases for RRab stars 
and the 1/2-day alias and first harmonic of the true period for 
RRc stars. 

This empirical overview of the most frequent spurious peri- 
ods, obtained by means of the simulated light curves, enabled us 
to quantify the relative importance of the different posible aliases 
present in our RRLS survey. Based on these results, an inspection 
of the identified period aliases was included in the period search, 
described in the previous subsection, during the visual inspection 
of the period-folded light curves. This allowed increasing the com- 
pleteness in the RRLS identification, which is described in the fol- 
lowing section. 

3.3.3 Final Selection of RRLS 

We visually inspected the light curves for the 66,086 stars for 
which a statistically significant (A ^ 3.0) minimum of Svri was 
found. The visual inspection and classification of light curves was 
done using the plotting tool InspectLC.py, developed as part of this 
work, whose graphical user interface is shown in Figure[8] The In- 
spectLC.py program allows visualizing, for each star, light curves 
in three different filters (in this case V,RJ) period-folded with each 
of the three most significant periods found. These light curves are 
shown in Fig. [founder the labels PI, P2 and P3. The text box to 
the left of the window summarizes each star's parameters: iden- 
tification number, a, 8, amplitude and number of observations in 
each of the three filters, the three best periods (Pj with i = 1,2,3) 
and their corresponding values of A. The 'TRY ALIAS' button dis- 
plays, for a given period Pj, the period-folded light curves using 
5 of the most common alias periods (see Table [5j plus one user- 
defined alias, computed from the values inserted in the 'k' and 
'1/p' textboxes (see Section 13.3.21 and Equation |5). These light 
curves are displayed in Figure [8] under the labels P4 to P8 and 
P9 respectively. The program does not display the light curves for 
alias periods outside the period range of RRLS. The 'PLOT P-A' 
button plots each period with its corresponding aliases in a period- 
amplitude plot, as shown in the upper right corner of the window in 
Figure [8] Thi s plot also shows data for RRLS in the globular clus- 
ter M3 from ICorwin & Carnevl J200lh as a reference to illustrate 
the loci ocuppied by RRab and RRc stars in this diagram. This 
plot is useful for deciding between different periods that produce 
similarly good period-folded light curves. In the above example PI 
was selected as the correct period. Aditionally, the 'DSS' button 
allows downloading and plotting a DSS thumbnail image centered 
on the coordinates of the star, in order to facilitate the identification 
of contaminants like spurious extended objects, artifacts on very 
bright stars or double stars. Finally, the 'SAVE' button allows sav- 
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Figure 7. Top: Recovered versus true period for synthetic stars: (Left) RRab (Right) RRc. Bottom: Loci defined by the period aliases visible in the top panels. 
The shaded red area indicates the identity locus. The remaining indicate the loci of the most common alias periods: (Green) one-day aliases, (Blue) 2 x P rea \ 
and 0.5 x P rea i, (Magenta) 2x one-day alias, (Cyan) 1/2 day aliases, (Orange) 1/3 day aliases, (Grey) 1/4 day aliases, (Yellow) 2x1 /2-day alias. The shaded 
areas have a width equal to 2 per cent separation from the corresponding locus. Note that, in the top panels, the points lying within the identity locus correspond 
to 80 per cent and 46 per cent of the total synthetic RRab and RRc samples respectively. 



ing for each star, in an ascii format file, upto two periods, comments 
on each period and the variable star classification selected from the 
bullets on the lower left part of the window. 

Finally, a total of 160 RRab and 51 RRc stars were identified. 
As a byproduct of the visual inspection, eclipsing binary candi- 
dates were also found and will be discussed in an upcoming pa- 
per. The catalogue of V,R and / light curves for all RRab and 
RRc stars is shown in Figure [A] in Appendix [A] The RRLS light 
curve parameters, such as amplitudes, maximum-light magnitudes 
and phase, were computed by fitting light c urve tem p lates t o the 
data using RRafc and RRc templates from iLavdenl i ll 9981) and 
IVivasetal.l ( l2008h respectively, and further refining the period with 
a step of 10~ 6 d. These parameters are listed in Tables [6] and [7] for 
RRa£> and RRc respectively. 

In the following sections we describe the effect of period 
aliases in the catalogue, the completeness estimated in the identifi- 
cation of RRLS, as well as the posible contaminants of the sample. 



3.3.4 Comparison with Variable Star Catalogues 

The prese nt survey partially overlaps with the RRLS catalogue 
from |V04| (see Figure [TJ, also conducted with QUEST observa- 
tions but focused on high latitude regions. All 21 RRab stars and 
31 (~ 78 per cent) of the RRc stars from V04, expected in our sur- 



vey's area were succesfully recovered. The remaining 22 per cent 
of RRc stars were re-classified as contaminants based on their im- 
proved light curves which include a larger number of observations, 
in particular in the R and / bands, with respect to V04. The 52 
RRLS recovered from V04 keep their original identification num- 
bers, which correspond to numbers ^ 999 in Tables[6]and|7] The 27 
RRLS with identi fication numbers in the range 500-600 are from 
IVivas et alj fepll in preparation). The 132 RRab and RRc stars 
identified in the present survey have identifiers from 1000 onwards. 

We also cros s-matched o ur RRLS catalogue wit h the ASAS-3 
jPoimanskil2002h and NSVS l lKinemuchi et all2006l) . with a toler- 
ance o f 7". Out of the 132 RRLS not present in V04 or lVivas et all 
42012b . 131 are new discoveries: one match is found in NSVS 
and none in ASAS-3. Our RRLS 1012 corresponds to star 722 
of the NSVS. The period reported for this star on the NSVS is 
Pnsvs = 0.53495 d which coincides within the errors with our pe- 
riod Pquest = 0.53497 d and in both catalogues, the star was clas- 
sified as RRab. 



3.3.5 Completeness 

The completeness achieved in the RRLS identification was esti- 
mated from the catalogue of synthetic light curves for RRab and 
RRc stars described in Sec. 13. II The completeness was computed, 
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Figure 8. Graphical user interface of the InspectLC.py tool developed for light curve visual inspection. 



Table 6. Light curve parameters for RRafo stars. 



ID 


a 


5 


(N v ^ R ,Nj) 




Amp V 


AmpR 


Amp/ 


HJD mnt7 


(V) 









(°) 


(°) 




(d) 


(mag) 


(mag) 


(mag) 


-2450000 (d) 


(mag) 


(mag.) 


(mag) 


1000 


60.935722 


-1.498512 


( 0. 10. 9) 


0.444400 




0.91 


0.62 


1200.70430400 




16.79 ±0.02 


16.55 ± 0.03 


hum 


62.136738 


-2.995186 


(0, 11,0) 


0.465461 




0.72 




2279.76109714 




19.16 ±0.11 




1 002 


62.444096 


-4.000918 


(27, 19, 11) 


0.498114 


1.16 


1.13 


0.87 


1823.74922626 


17.53 ±0.04 


17.13 ±0.02 


16.49 ± 0.02 


1003 


65.498161 


-3.226088 


( 28, 16, 24) 


0.501271 


1.25 


1.07 


0.87 


1824.19917603 


15.91 ±0.01 


15.77 ±0.02 


15.40 ± 0.01 


563 


66.183346 


-3.859660 


(24, 0, 11) 


0.539828 


0.85 




0.41 


1867.80606764 


14.03 ± 0.02 




13.60 ±0.01 



Note: Table[6]is published in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its form and content. 



independently for each RRLS type, as the percentage of synthetic 
stars recovered after the period search described in Sec. 13.3. ll was 
performed on a sample of ~ 65,000 synthetic RRLS of each type. 
A synthetic star was considered succesfully recovered if one of the 
3 best periods found were within 10 per cent of the true light curve 
period or one of the period aliases described in Sec. 13.3.21 The av- 
erage completeness was computed for synthetic stars in the magni- 
tude range 14.0 ^ V ^ 18.5 as a function of position in the survey 
area, since it strongly depends on the time sampling and number of 
observations in each filter. 



The resulting completeness maps for RRLS of types ab (left) 
and c (right) are shown in Figure [9] in equatorial coordinates. As 
illustrated by the figure, for RRab stars the average completeness is 
~ 95 per cent for most of the survey area, except in the declination 
stripe —6° ^ 8 ^ —4°, where the number of observations per object 
is typically less than 15 in all three photometric bands (see Figure 
[5}- For RRc stars the estimated completeness is ~ 80 — 85 per cent, 
significantly lower than for RRab stars, as expected since RRc light 
curves have smaller amplitudes. 



3.3.6 Possible Contaminants 

The main source of contamination in an RRLS survey are eclips- 
ing binaries (mostly of the WUMa type), 8 Scuti and SX Phoeni- 
cis stars. Eclipsing contact binaries or WUMa stars, consist of 
a system in which both stars have filled their Roche lobes and 
are undergoing mass transfer dSterken & Jaschekll2005r) and have 
light curves showing two eclipses of very similar depths. 8 Scuti 
stars are pulsating main sequence dwarfs and blue stragglers, 
while SX Phoenicis stars are their metal -poor counterparts; both 
types are also called Dwa rf Cepheids dNemec&Mateolll99rj 
ISandage & Tammannl l2006h and have asymmetric light curves, 
similar to those of RRLSs, though w ith smaller periods (P < 0.2 
d) and amplitudes (AmpV < 0.2 mag. lSterken & Jaschek 2005). 

A plot of V amplitude versus period is shown in Figure [T0l 
for RRab, RRc, Algol, [3 Lyrae, W UMa and 8 Scuti stars from the 
ASAS-3 survey l lPoimans M2002I) . As illustrated by the figure, the 
period and amplitude ranges for most of these stars largely overlap, 
making the light curve shape an important criteria in distinguishing 
between the different types of variables. Most of the contamination 
problems affect the RRc rather than the RRab stars. 

Eclipsing binary light curves, of any type, are unlikely to be 
mistaken for RRLSs of type ab, while dwarf cepheids, though hav- 
ing similar light curves, have periods and amplitudes which are 
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Table 7. Light curve parameters for RRc stars. 



ID 


a 


5 


(N v ^ R ,Nj) 




AmpV 


AmpR 


Amp/ 


H1D maxl 


(Y) 


<«> 







(°) 


(°) 




(d) 


(mag) 


(mag) 


(mag) 


-2450000 (d) 


(magi 


(magi 


(magi 


1013 


74.151787 


-3.740614 


(21. 17. 11) 


0.331910 


0.43 


0.40 


0.23 


1824.12209720 


17.91 ±0.05 


17.67 ± 0.04 


17.46 ± 0.05 


1015 


74.725197 


3.882039 


( 44, 39, 33) 


0.400277 


0.58 


0.36 


0.23 


1870.98566120 


15.22 ± 0.01 


14.83 ± 0.01 


14.48 ± 0.01 


1 1 '■> 


75.680612 


5.070141 


( 19, 12, 19) 


0.384742 


0.38 


0.38 


0.33 


1871.09802948 


15.90 ±0.01 


15.77 ±0.01 


15.54 ± 0.03 


1022 


76.733847 


0.805251 


( 33, 20, 9) 


0.412460 


0.44 


0.37 


(1.16 


1539.59914360 


18.80 ±0.07 


18.63 ± 0.06 


18.38 ±0.10 


1028 


77.485196 


0.078322 


(82,21,34) 


0.370848 


0.30 


0.12 


0.14 


1539.62626816 


17.51 ±0.04 


17.40 ± 0.03 


17.21 ±0.05 



Note: Table[7Jis published in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its form and content. 
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Figure 9. Map of average completeness in the identification of RRafo (left) and RRc (right) stars, computed from synthetic RRLS within the photometric 
completeness limits of the survey (14.0 ^ V ^ 18.5). 



too small for RRab stars. A few dwarf cepheids could have alias 
periods in the range of RRLS, however in this case the contami- 
nants are unlikely to follow th e well known Period- Amplitude re - 
lationship for RRab stars (e.g. ISmith|[l995l : ICacciari et alj|2005h . 
which can in turn be used as a criterion to discard them. This cri- 
terion is also helpful to avoid other possible contaminants of the 
RRab sampl e, such as Anomalous Cephe ids and short period Type- 
II Cepheids dSandage & Tammannll2006h . also known as BL Her- 
culis and W Virginis stars, which have light curve shapes similar 
to RRab stars with periods in the range from 1 — 50 d (W Vir) and 
1 — 7 d (BL Her), partially overlapping the range for RRab stars in 
the long-period end. Therefore, the contamination expected in the 
RRab sample can be considered negligible. 

The RRc sample is, on the contrary, more prone to con- 
tamination. The light curves of eclipsing binaries of the WUMa 
and p Lyrae type can resemble those of RRc stars if period- 
fol ded with half the true per iod of the binary. However, as noted 
by iKinman & Brownl ( feOlOh . an important difference between 
RRc stars and eclipsing binaries is that pulsating stars have larger 
amplitudes in bluer bands compared to red bands; while eclips- 
ing binaries have similar amplitudes in different ph otometric bands 
dKinman & Brownll201filSterken & Jaschekll2005h . This is a use- 
ful criterion for differentiating both types of variables, except for 
variable stars with poorly sampled light curves or very small am- 
plitudes (< 0.2 — 0.3 mag), in which case the distinction is more 



difficult. On the other hand, since this behaviour of the amplitude 
in different photometric bands is common to pulsating stars, it is 
not useful for distinguishing RRc stars from the larger amplitude 
8 Scuti or SX Phoenicis stars, making these likely contaminants of 
the sample. 

Therefore, the contamination expected in the RRc sample 
could be significant. In particular, among the RRc, a large number 
of WUMa and S Scuti stars belonging to the thin disk is expected 
given the low latitudes being surveyed and the fact that we are not 
impossing any colour cuts which wou ld sizeably have reduced the 
contamination dlvezic et alj200d.lv04l) . This added to the relatively 
low completeness of the RRc sample (~ 80 per cent), makes it in- 
appropriate to use the RRc stars in the computation of the thick disk 
structural parameters, which is why we have restricted the follow- 
ing analysis to the RRab sample. 



3.4 Physical Parameters of RRL stars 

3. 4. 1 Extinctions 

The present survey spans areas at very low galactic latitudes (|£>| < 
15° ), where the extinction can be very h igh and spatially variable. 
The lSchlegel. Fin kbeiner & Davis! dl998l) dust maps, although cov- 
ering the whole sky, are unreliable at galactic latitudes lower than 
~ 10°. On the other hand, a correct estimate of the extinction is 
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Figure 10. Bottom: V amplitude versus period for variable stars of type 
RRqfc RRc, 8 Scuti, A lgol, |3 Lyrae and WUMa from the ASAS-3 cata- 
logue (Pojmanski 2002). Top: Period histogram for the different variable 
star types illustrated in the bottom panel. 



crucial in the determination of distances for the survey RRLSs and 
their subsequent use in the computation of density profiles. 

In addition to being excellent standar d candl es, RRLSs can 
also be used as colour-standards. ISturchl dl966h finds the av- 
erage B — V colour index during the phase of minimum light 
((]) £ [0.5,0.8]) is approximately constant for RRa£> stars, having 
a small dispersion of ~ 0.03 mag; which makes it possible to 
estimate the extinction affecting a particular RRab star by sim- 
ply comparing the observed minimum-light colour with the in- 
trinsic colour. L ater on, the calibration was ext ended to V — I by 
Day et al.l j2002h and IGuldenschuh et al.l d2005l) . and to V-R by 



Kunder. Chabover & Lavde d fcOld) ; obtaining smaller standar de- 



viations of ~ 0.02 mag in these colour indices. 

In the present work, we used the mini mum-light colour cal- 

0.5-0.8 



ibrations (V — R /m 

and(K-C r °- 8 = 
For the 

we used the intrinsic colour (R — I)^„ 7 ' 8 



0.27 ± 0.02 from iKunder et al.l fefJlCh 
0.58 ±0.02 from IGuldenschuh et al.l j2005h . 
40 RRab stars that do not have V-band observations 

0.31 ±0.03, com- 
puted from the previous two colours. The observed minimum-light 
colours were computed from the subtraction of the best fitting 
templates in the two bands used. Comparing these with the as- 
sumed intrinsic colours we computed the corresponding E(V— R), 
E(V — I) and E(R — I) colour excesses and extinctions Av,Ar 
y A j were derived assuming the standard reddenning law from 
ICardelli. Clayton & Mathisl dl988h . 

We used the 56 stars with observations in all three filters, to 
compute the residuals in the E(B — V) derived using the V — R, 
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Figure 11. Residu als between colour excesses E(B — V) from 
ISchlegel et ail < 19981) dust maps and those computed from RRab in- 
trinsic colours, as a function of galactic latitude b. 



V—I and R — I colour indices, in order to test the consistency of the 
resulting extinctions. The mean differences obtained were E (B — 
V ) VR - E {B - V ) vi = -0.007 mag and E (B - V) m - E {B - V ) y/ = 
—0.006 mag, with standard deviations of 0.092 mag and 0.082 mag 
respectively, showing there is very good agreement between the 
different measurements. For stars with data in the VRI filters the 
adopted E(B — V) were computed as the error- weighted mean of 
the colour excesses computed from the K — R and V — I intrinsic 
colours and the corresponding errors were computed by standard 
error propagation. 

Residuals between the se colour excesses an d those obtained 
from interpolations in the ISchlegel et al.l q 1998T) dust maps are 
shown in Figure QT] The plot shows the residuals are nearly zero 
for moderately high galactic latitudes (\b\ ^ 20°), while as the lat- 
itude decreases the colour excesses computed from RRafc intrinsic 
colour s are systematically smaller than those from lSchlegel et all 
dl998h . This behaviour is expected since the former correspond to 
the colour excesses integrated up to the distance of each RRLS, 
while the latter are integrated along the entire line of sight. Fig- 
ure QT] also shows larger residuals around latitudes b ~ —20° 
and b ~ ±10°, which are presumably due to abrupt variations in 
the extinction on relatively small scales, mainly due to the pres- 
ence of the Orion Molecular Cloud (at b ~ —20°). The points 
showing the larger residuals in Figure QT] coincide with regions 
where the extinction changes abruptly and the colour excesses 
derived f rom RRab star s are s ystematically larger than those re- 
ported bylSchlegel et al.l dl998h . which is to be expected since the 
ISchlegel et alj maps have relatively large pixels with 6.1' x 6.1' 
and abrupt variations on smaller scales will be averaged out, while 
RRaii stars probe the extinction exactly on the line of sight towards 
each star. 



3.4.2 Photometric Metallicities 

RRLSs of type ab exhibit a well known relationship between metal- 
licity [Fe/H], pulsation period P and t he (foi phase of a Fourie r 
light curve decomposi tion, discovered bv ljurcsik & Kovacsldl996l) . 
Also, ISandagel d2004|) finds that <|>3i is correlated with the light 
curve amplitude and fi nds an expression analogous to that of 
ljurcsik & Kovacj d 19961) relating [Fe/H] with the light curve am- 
plitude and period for RRab stars. 
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Figure 12. £"(B — V) map from lSchlegel et aljjl998l) . Open circles indicate survey RRab stars. 



3.4.2.1 Jurcsik & Kovacs Metallicities ljurcsik & Kovacsl 
dl996h found an empiric relationship between the metallicity 
[Fe/H], period P and <|>3 1 phase of a Fourier decomposition of 
RRab light curves, expressed in their Equation 3 as 



[Fe/H] = -5.038 -5.394P + 1. 345031 



(6) 



The dispers i on in this relationship was estimated by 
ljurcsik & Kovacsl ( 1 19961) to be 0.11 dex from a comparison with 
spectroscopically determined metallicities. 

The Fourier fitting process is very sensitive to the light curve 
sampling. For relatively well sampled curves, but with a moder- 
ate number of observations (N ~ 30), the fits usually have exces- 
sive ripples which translate into undesirable uncertainties in the 
determination of (j>3 1 . In ord er to circumven t this problem, we 
used the method proposed by iKovacs & Kupl d2007l) which con- 
sists in fitting an observed light curve with one of their 248 light 
curve templates, representative of RRab stars, and determining 
(j>3 1 from the Fourier decomposition of the best-fitting template 
This allows for robust Fourier decompositions of light curves with 
a moderate number of observations (N > 20). In order to have 
an estim ate of the quality of the fit, we use the Df index de- 
fined by ljurcsik & Kovacsl l ll996l) . which tests for consistency of 



6 We used the tff.f software from [Kovacs & Kupi 1 2007b in the template 
fitting and Fourier decomposition processes 



the Fourier decomposition parameters using the fact that different 
Fourier phases (|)y are correlated (see their Table 6). The restric- 
tion Df < 3 was applied in order to ensure the reliability of the 
Fo urier fits used in the com putation of metallicities, as suggested 
bv ljurcsik & Kovacsl dl996h . 

3.4.2.2 Sandage Metallicities Using a samp le of RRab stars 
having spectroscopically derived metallicities, ISandagel d2004l) 
finds the following empirical relationship between [Fe/H] and the 
light curve period, amplitude 



[Fe/H] v = (-1.45 ± 0.04) Amp v 



(7.99 ±0.09) log P 

-(2. 145 ±0.025) (7) 



ISandagel argues this relationsh ip reflects the same behaviour 
found by Jurcsik & Kovacsl dl996l) . since (foi corre lates directly 
with the light curve amplitude. In practice. ISandager s relation can 
be used in many more RRLS catalogues since it only requires an 
accurate knowledge of period and amplitude, w hile in order to per- 
form the Fourier fits o fljurcsik & Kovacsl d 19961) the full time series 
data are required and the light curves must have an appropriate time 
sampling. 

The relationship expressed in Equation [7] is calibrated for the 
V band, found in the majority of RRLS observations in the litera- 
ture. In our catalogue of 160 RRab stars, 25 per cent of the stars 
have observations only in the R and/or / bands and not in V, which 
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Figure 13. Distribution of differences for photometric metallities [Fe/H]/; 
and [Fe/H]/ computed from R and / band amplitudes (Equation \§\ and |5J, 
with respect to [Fe/H]y, computed from the V band amplitude (Equation 



made it necessary to extend Sandage's relation for these photo- 
metric bands. Using the photometric metallicities computed from 
Equation |7] for 49 of our stars which are well-observed in V and 
R (AV > 20 and N R > 20) and 45 stars in V and / (N v ^ 20 and 
Nj ^ 20), we obtained analogous equations for the R and / bands, 
via least squares fitting of the corresponding metallicities, periods 
and amplitudes. The corresponding equations are 

[Fe/H]/} = (- 1.052 ± 0.027) Amp/} 

- (6.281 ±0.087) log P- (2.254 ±0.028) (8) 

[Fe/H]/ = (-1.469 ±0.032) Amp, 

- (6.516±0.091)log/>- (2.250±0.027) (9) 

Figure Q~3] shows the distribution of residuals of photomet- 
ric metallicities [Fe/H]/} and [Fe/H]/ computed from the R and / 
band amplitudes and Equation[8]and |9]respectivel y, with respect to 
[Fe/H]y derived from Equation[JJ(Sandage 2004). The two distri- 
butions show the agreement between the photometric metallicities 
computed from Equation [8] an d [9] derived in th is work, and those 
computed via Equation 171 from [Sandagej d2004l) . resulting in mean 
residuals of —0.01 dex and 0.006 dex respectively. The intrinsic 
dispersion of photometric metallicities computed from Equation 
|7j in comparison w ith spectroscopic metallicities, was estimated 
by ISandagel d2004l) as c |p e y H ] = 0.26. Standard deviations for the 
residual distributions in Figure Q~3] were found to be 0.19 dex and 
0.17 dex for metallicities computed from R and / band amplitudes 
respectively. Adding these dispersions in quadrature with Ofc e/ < H ] , 
we estimated the intrinsic dispersion of Equations [8] and |9j to be 
rjpe/H] = 0.32 dex and CT [p e /H] = 0-31 dex respectively. Uncertain- 
ties in the reported photometric metallicities were computed by 
adding in quadrature the errors obtained from error propagation in 
Equations |7|9| with the respective intrinsic dispersion of each equa- 
tion. Using these procedures, the first term accounts for the ampli- 
tude uncertainty (the contribution of period uncertainties are negli- 
gible) while the second term accounts for the intrinsic dispersion of 
the method. The typical uncertainties obtained for the photometric 
metallicities are ~ 0.4 dex. 
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Figure 14. Distribution of residuals of photometri c [Fe/H] f ol and sp ectro- 
scopic [Fe/H]™; metallicities for 27 RRab in the I Vivas et al J JSoOSh cata- 
logue. 



An additional test on these methods was performed by com- 
paring the s pectro scopic metallicities for 27 RRab stars from the 
IVivas e t al. (2008) catalogue, with the corresponding photometric 
metallicities computed from Equation [7J Figure [14] shows the re- 
sulting distribution of residuals [Fe/H]y of — \Fe/H] spc . The mean 
and standard deviation of these residuals are —0.01 dex and 0.34 
dex respectively, in agreement with the intrinsic dispersions for 
Equation[7J[8]and[9] 

3.4.2.3 Metallicity Distribution The photometric metallicities 
adopted for RRab stars were calculated from a weighted average 
of the metallicities computed using the ljurcsik & Kovacsl J 1996b 
(Equatioi© and lSandageH2004l) (Equation l7l9l > methods. For those 
stars for which the Fourier decomposition obtained w as unreliable 
(i.e. Dp > 3), the metallicity computed from lSandagej ' s method was 
adopte d. The metallicities reported here are on the IZinn & West! 
i 19841) scale, and for all stars lie inside the metallicity range 
(— 2.5 [Fe/H] ^ ±0.07) of th e stars used in the calibration of 
the ljurcsik & Kovacsl i 19961) and lSandagei J2004I) methods. 

The distribution of metallities obtained is shown in Figure 
|15| for the full RRab sample (black) and for the stars close to the 
Galactic Plane with |z| ^ 2.5 kpc. The total metallicity distribution 
extends from [Fe/H] ~ —2.5 up to [Fe/H] ~ ±0.5 with a mode 
at [Fe/H] ~ —1.48 as can be seen in the figure. The main peak 
coincides with the expected peak of the halo metallicity distribu- 
ti on, which has an obs erve d metallicity [F e/H] = — 1.5 according 
to lCarnevetaTl i ll 9941) and lLavderJ l ll995h , for F and G field halo 
dwarfs and RRLs respectively. The metallicity distribution near the 
Galactic Plane shows a hint of an excess for [Fe/H] ^ — 1.0, around 
the metallicity expected for thick d isk stars around [Fe/H] ~ —0.8 
dGilmore et al]l995l ; lLa"vderJI 19951) . 

3.4.3 Distances 

The computation of heliocentric distances for the RRab stars was 
done using absolute magnitudes My RL an d M? RL derived from the 
Period -Luminosity-Metallicity relations in lCatelan. Pritzl & Smith! 
J2004l their Equat ion 8 and 3 respectively). According to 
Cate lan et alj <2004h the use of Period-Luminosity-Metallicity re- 
lations is preferable over My — [Fe/H] relations, since the former 
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Figure 15. Photometric metallicity distribution for the 160 RRafo stars in 
the survey (black). The filled histogram (blue) shows the distribution of 
metallicities for stars near the Galactic Plane with \z\ ^ 2.5 kpc. 
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4 SUMMARY 

We have presented the QUEST RR Lyrae survey at low galactic lat- 
itude, which spans an area of 476 deg 2 on the sky, with multi-epoch 
observations in the V, R, and / photometric bands for 6.5 x 10 6 stars 
in the galactic latitude range —30° ^ b ^ +25°. The survey has a 
mean number of 30 observations per object in each of the V and 
/ bands and ~ 25 in the R band, while in the best sampled areas 
each object can have up to ~ 120 — 150 epochs in V and I and 
up to ~ 100 in R. The completeness magnitudes of the survey are 

V = 18.5 mag, R = 18.5 mag and / = 18.0 mag with saturations of 

V = R = 14.0 mag and / = 13.5 mag. 

The survey has identified 211 RRLSs, 160 bona fide stars of 
type ab and 51 strong candidates of type c. The RRLS catalogue 
presented here contains the positions, mean magnitudes in VRI, pe- 
riods, amplitudes, photometric metallicities, distances and individ- 
ual extinctions computed from minimum light colours for each star, 
as well as the period-folded light curves for RRab stars (Tables [6] 
[7] [8] and Figure [A}. The typical distance errors obtained were ~ 7 
per cent. The completeness of the RRLS survey was estimated in 
> 95 per cent for RRab and ~ 85 per cent for RRc stars, and was 
computed from light curve simulations reproducing the time sam- 
pling and other observing characteristics of our survey, as well as 
the variable identification and period searching techniques used. 

Our RRLS survey spans simultaneously a large range of helio- 
centric distances 0.5 ^ i?/ lf /(kpc) ^ 40 and heights above the plane 
— 15 ^ z(kpc) ^ +20, with well characterized completeness across 
the survey are a. Combining ours with a bri ght RRLS survey (such 
as those from iLavdenlll 9951 : iMain tz 2005, see Sec.QJ to increase 
the distance coverage even further, will result in a very good tool for 
studying the Galactic thick disk's structure, in particular for the de- 
termination of the density profile's scale length, scale height as well 
as their dependence with radial distance in order to explore the flare 
and the truncation profile of the thick disk. This way, the QUEST 
RRLS from this work will help put constraints on the structure of 
the more external parts of the thick disk toward the Anticenter, as 
well as the mechanisms that could have contributed to the forma- 
tion of the Galactic thick disk. These issues will be addressed in 
detail in an upcoming paper of the series. 



Figure 16. Distribution of heliocentric distances of survey RRofe stars. 



allow accounting for evolution off the Zero- Age Horizontal Branch 
(ZAHB). The total m etal abundance Z was computed from [Fe/H] 
using Equation 10 in ICatelan et all j2004h . assuming an a to Fe 
abun dance of [rx/Fe]=0.2 dex typical of h alo and thick disk stars 
(e.g. I Venn et al.l [20041 ; iReddv et alj|2003l) . For the 9 RRab stars 
which only have R band photometry, we computed M§ RL by inter- 
polatin g the corresponding coefficients in Table 4 of ICatelan et al.l 
( 120041) . 

Finally, the heliocentric distances were computed from Pog- 
son's equation log/J^; = (m,- — A; — M, +5)/5, where A/, M; and m, 
correspond to the extinction, observed and absolute magnitudes in 
the i filter. The corresponding errors were computed via error prop- 
agation in this equation, resulting in typical distance errors of ~ 7 
per cent. The resulting heliocentric distance distribution is shown 
in Fi sure [76l and spans the range 0.2 ^ i?/, e /(kpc) ^ 80. In Table[8] 
we summarize the colour excesses, metallicities, heliocentric dis- 
tances Ri le i, as well as height above the Plane z and radial distance 
R projected on the Galactic Plane, obtained for all RRab stars. 
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Table 8. Distance, colour excess and metallicities for KRab stars. 



ID 


E(B-V) 


[Fe/H] 


Rhel 


z 


R 




(mag) 


(dex) 


(kpc) 


(kpc) 


(kpc) 


1000 


0.02 


-1.0 ±0.3 


17.1 ± 1.3 


-10.4 ± 0.8 


21.5 ± 1.0 


1001 


0.13 


-0.9 ± 0.3 


46.4 ± 2.5 


-28.0 ± 1.5 


44.8 ± 1.9 


1002 


0.41 


-1.4 ±0.5 


13.0 ± 1.7 


-7.9 ± 1.0 


18.1 ± 1.3 


1003 


0.01 


-1.6 ±0.1 


11.5 ±0.8 


-6.5 ± 0.5 


17.3 ±0.6 


563 


0.01 


-2.1 ±0.1 


5.0 ±0.2 


-2.8 ±0.1 


12.0 ±0.2 



Note: Table[8]is published in its entirety in the electronic edition of the journal. 
A portion is shown here for guidance regarding its form and content. 
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APPENDIX A: RRLS LIGHT CURVES 

The catalogue of V, ^ and / period-folded light curves for all 
RRab and RRc stars found by the survey, is shown in FigurelAl. 
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Figure Al - continued Light curves of survey KRab and RRc stars. 
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Figure Al - continued Light curves of survey KRab and RRc stars. 
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Figure Al - continued Light curves of survey KRab and RRc stars. 
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Figure Al - continued Light curves of survey KRab and RRc stars. 
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Figure Al 



continued Light curves of survey KRab and RRc stars. 
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Figure Al - continued Light curves of survey KRab and RRc stars. 



